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ABSTRACT 



Relationships relating power, magnetic field, current density, and 
ampere-turns in terms of certain dimensionless factors are summarized 
for many types of coil geometries and current distributions. A num¬ 
ber of plots of these factors are presented. 

The field homogeneity in magnet structures is presented in terms 
of a series expansion about the origin utilizing Legendre Polynomials. 

A number of tables to facilitate design of homogeneous fields are 
presented. A method of achieving homogeneity in long solenoid struc¬ 
tures by the use of determinants is discussed. Expressions for the 
axial field from uniform and radially varying current density coils are 
given. 
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PART I 

RELATIONSHIPS BETWEEN MAGNETIC FIELD, POWER, AMPERE - TURNS 
AND CURRENT DENSITY 

It has been customary In the literature on magnet design to write the relationships 
between fields, power, ampere-turns and current density In terms of a number oi dimen¬ 
sionless factors. In this section we will treat In summary form twelve cases of coll 
construction (Cases I - XII). This work is In part a summary, and enlargement of parts 
of the following three papers: F. Bitter, R. S. L Vol. 7, 1936 Part II; W. F. Gauster, 
AIEE, Fall General Meeting, Chicago, October 1959; F. Gaume, Journal de Recherche 
du Centre Nationale de Recherche Scientifique, No. 43, June 1958. (English translation, 
Lincoln Laboratory, M81-12, part I and II, Dr. H. H. Kolm.) 

A. The most commonly used dimensionless factor Is the G factor first suggested 
by Fabry, which connects a given magnetic field with the power required to produce this 
field and la a function of the type of current distribution and the geometry of the coll. 

-°(sr 

with W watts 

p - resistivity In ohm - cm 

Ql = inner radius In cm 

X - fractional volume of conductor 

The space factor A la assumed constant throughout the volume. If A is a func¬ 
tion of x and y, t>ore complicated formulas are needed. 

B. A seco v dimensionless factor J has been suggested by W. F. Gauster and 
connects the current density at the innermost winding of a magnet to the power of the 
magnet. J Is a,V-.. . function of the current distribution and the geometry of the coil. 



C. i v. ..itlcular type of coll a maximum G factor can be found which will 

give the most magnetic field for the least power used. However, in many types of coll 
construction this maximizing of the G factor may lead to excessively high current densi¬ 
ties at the innermost winding which may make cooling the coil difficult. To illustrate 
this point a ratio, y , can be formed which compares the maximum current density In 
any coil to the current density In a uniform current density coil, where both havo the 
dame Inner radius and both develop ihe same field. Quantities for the uniform case are 
designated as primed: 







(3) 


if H=H' 

/wf G' 

[WJ “ ~G (3 > 

from j/wV* 

y x ~ r U 7 

and using (3) 

j,' ~ r g y (4) 

This ratio indicates the difficulty of cooling and supporting any coil compared with 
a coil of uniform current density at the same field and bore. 


D. The ratio in equation <3) is also useful and represents the power saving in 
using any current distribution compared with uniform current density when generating 
the same field in the same bore. We give it in the more useful form. 



These dimensionless factors and ratios are presented in Table I for nine types of 
magnet construction aud following Table I for three additional more complicated magnet 
constructions. 

One of the most Interesting new coil constructions presented is that of case V. 

This coll was first suggested by F. Gaume and achieves a relatively high G factor with¬ 
out excessive current densities and can be quite simply constructed of plates In the 
usual Bitter magnet style, with the introduction of axial current variations by mak ing the 
plates thicker towards the ends of the magnet. 

A number of plots of the geometry factor, G, against CT and 0 axe given In 
Figures 1 through 6. Cockcroft's curves for case VIII are reprinted in Figure 1, Bitter's 
curves for case VI in Figure 2, new curves for cases VIII and VI in the small CT and 
0 region are given in Figures 3 and 4. New curves for case V (Gaume's axial current 
variation) and for case X (the ellipsoidal shell) are given in Figures 5 and 6. The 
variables <X and 0 are defined: 











_L 


-Eb- 


1 

§§l§g 




a- a */a, 


% 


at 


E. It Is sometimes Instructive to write the magnetic field In terms of ampere- 
turns rather than power. This irises from the practical standpoint that It is difficult 
to always predict the resistance of a magnet exactly, especially as a function of power, 
and one may wish to know how much field will be produced for each ampere available In 
the power supply. This relationship between ampere-turns and field can be related to 
the J factor and the G factor an! miscellaneous CC and /3 terms and numerical con¬ 
stants. The relationships are given here for three cases: Case VI, Case VII and Case 
VTH. All the others are easily derivable by integrating the current density over the coil 
to obtain the total current in terms of J and W. Total current I and the factor J can now 
be substituted for W in equation (1), and the ampere-turns related to the field. 


Case VI: 


Case VII: 


Case VIII: 


H 6 


IA[_ _L_] 

t UiAncrJ 


G 6 
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( 8 ) 


?. Useful geometry factors can be written for very long solenoids (long enough so 
that end effects can be neglected). 

For a constant current density coll the ampere-turnB per unit length can be 
written: 
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NI ^ /(W/2b) A \ 1/2 / lMV /2 
2b “ \ 2^yo / \ a* 4-1 / 

where W/2b is the power per unit length (watts/centimeter) 

For the disk construction, with a radial current distribution we write 


fSW'M' 


Rewriting (9) and (10) in terms of the field in gauss and power per length as watts 
per meter, we obtain f 


H (gauss) = G s ( ^tb/raekrl x y^ 

50 \ p (ohm cm.) / 

(3 5 (uniform) = . fel l^ 


G 5 (radial) = (im cr) 1 ' 


These factors are plotted In Figure 7 and it is interesting to note that below an CT 
of 2. 5 there is negligible power saving with the disks, but at greater <2" , the saving is 
increasingly pronounced; at an CT of 10, 45% more power would be required by the uni¬ 
form wound solenoid for the same magnetic field. 

G. It Is useful to use the G factor notation when deallig with split coils (colls 
separated by a gap). If we take a pair of coils with parameters as defined in Figure 8, 
we can write the field In the gap by simple addition and subtraction of coils. 
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Figure 8 


(ii) 





















CASE X (FIGURE 6) 

Coll Shape: Ellipsoidal shell, CC = constant 



Geometry Factor: 


i) *>« 

' J 3 > 
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X3<1 

cr>0 
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/3>f 
<r >/3 


_J__ f (/SH/S'-cr*)'*)) 

( 0 Z -^W U ( ~S )) 

r ~ ' Ir ' ;Z n X M J_ 1 1 t 

5/3 \(/ 3 z -i ) V2 sec ii ~{a 2 -j 3 2 ) vz ** 


1 _ 

(a 2 -/ 3 2 )‘ /z 



G mm = .204 at a- = 6 ,/ 3 -Z 


Current Density Factor J: 


/—i_y* 

\ 4 7 r/i Xn < 2 * J 


i 4 7r/i in c 
at cr-6,/3=~~2, G^x ~ 0.204 
and Jfncu ** 0.149 


Current Density Itatlo ~y for samo Field referred to case VIII: - 1.31 
Powei Itatlo P for same Field referred to case VIII - 0. 76S 




CASE XI 


Coll Shape: Cylindrical Bore, tapered ends, non-coinciding apexes. 



b = ky-f-o 
k = slope 



X 


Geometry Factor: 


r ^ -7T</g _1_f k 

5 k(cr-0+7^10’ Lk 2 +1 

n ((k"+l , )a’ 2 +2kycr + y 2 ) , / 2 + a’ (k 2 +t) + y(b^k z +1 

((k 2 +1) + 2ky+y 2 )>/2 4^ 

+ ln fcr (Ck z +0-b2ky-ky 2 ) >/2 + y + k _ 

\ (Ck 2 f l) a 2 -f 2kycr +y 2 ) ‘/2 + y + W/J 


G is a maximum when k = 1 and CC =4.5 which reduces case XI to case VII Current 
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Density Factor J: 








CASE xn 

Coll Shape: Cylindrical Bore, tapered ends, non-colnclding apexes. 



b - ky+-c 
k * slope 


0 2 




tr** k VI _ 1 \ 

5 (W)^LVk(a*-l) + 7Wm_\/ # 

\y+<rkr 

i/z f VA 

^ ((k z +l)a- 2 + 2kya-fy 2 )+cr(k z +0 + y\l^TT/ 1 

((k 2 +t) +2ky + y ^'^+(F + i) ,/2 + r 


G Is a maximum when k = 1.0 and <X = 4.5 which reduces caso XII to case VII Current 
Density Factor J: 

1 , \'/z 


4ir + y 







L/R = a, \/p. 10.B=G,y^or 

Reproduction of Cockcroft’s G-factor curves for uniformly 
wound coils with a square cross-section (Case VIII) 
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A. 


PART H 

HOMOGENEOUS MAGNETIC FIELDS 


A magnetic field cm be written as a power series involving I^gendre polynomials 
expanded about the origin of the field. If this origin is on the axis and on the plane of 
symmetry the expansion will have no odd terms and the axial component Hz and the 
radial component Hyo can be written: 

H * (t '^ “J,T2 hW r2 " 2 P 2n-2 W 


U- COS & 


In these equations the coefficients are defined as follows: 

,.(2n-2), . d 2nZ H z (z,o) 

Hz (z ‘°' = — - dz 2n-2 - 

dz Jz = 0 

and P^n-2 (u) are tlie Legendre polynomials tabulated In Table 2. 


The Importance of this approach, due to W. Garrett 1 , Is that the field at any point 
on the axis of a system of coaxial coils can be given explicitly In simple algebraic form. 
From this., the derivatives of equations (1) and (2) can bo obtained, and so a complete 
expression for the field In the vicinity of the origin can be obtained. If sufficient deriv¬ 
ative terms are used the field can be found anywhere within a sphere whose center lies 
at the origin and which does not extend far enough to Include any comers of the coil. How¬ 
ever, out to within a few percent of the inner radius of the coil, the field can be found 
quite accurately with only a few derivative terms. Compensation can be designed for each 
component of the field separately. In general. If we have a coll system of the assumed 
symmetry having N variable parameters such as size, or ampere-turns, it is possible to 
specify N coefficients in the expansions 1 and 2. 


For example, to examine the standard Helmholtz pair of coils, where the colls are 
so spaced that there is no second derivative of the field at the origin, we can think of 
the coil system In the following way: 


1 . 


M. W. Garrett, J. Appl. Phys. 22, 9, Sept. 1951 
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a » d “G 

Figure 1 

As shown in Fig. 1, we consider a single coll having length equal to the coils and the 
gap separation, from which we subtract a second coil of the length of the gap. The gap 
coil is thought of as having a current in the reverse direction to the original long coil. 
We now merely need select the length of the gap coll so that it will have a second deriv¬ 
ative equal and opposite to that of the original long coll and the two will thus cancel at 
the origin. 

There are many other combinations which would also cancel the second derivative 
at the origin, for Instance, the hypothetical gap coil with its reverse current could have 
been somewhat longer and had less reverse current in it, or shorter and had more re¬ 
verse current In it. Both of these solutions, of course, lead to a continuous solenoid 
with no gap but with an altered current density in the center section. These various 
solutions will have both an effect on the amount of field in the center of the magnet and 
on the size of the higher derivatives. The greater the distance over which one compen¬ 
sates, the more parallel will be the flux lines from the compensating coil, hence the 
flatter will be the flux distribution at the origin, and hence the smaller will be the sum of 
the higher derivatives. This compensation may not be compatible with obtaining the 
highest central field for the least power, howeverf 


Equations (1) and (2) can be rewritten In the following formr^ 

H z (r,o) - H z (o,o) [h -£ 2 (^'>2 {>*) + <£* (|j) P 4 (u) + —] < 3 > 

H^Cno)- h z (o,o) [o+£ z ^) 2 Pi ( u )+£ 4 (a;) p 4( u )+•"] (4) 

°(sr 

2. R. S. Ingarden and J. Michalczyk, Bulletin De L’Academie Polonaise Des Sciences 
(Serle des sci. math, aatr., and phys) VIU, 5, 1960. 






The expansion has been normalized to the Inside radius of the coils a^. A tabulation of 
the coefficients for three useful geometries is given in Table 1. The geometries 
treated are (i) a cylindrical thin current sheet, (li) a finite solenoid with a uniform 
current density, and (iii) a finite solenoid with a radial current distribution )• 

The value of the first eight Legendre polynomials and their derivatives are given for con¬ 
venience in Table 2. From equations (3) and (4) the magnetic field can be found at any 
point near the origin of a magnet, at a distance T from the origin and at an angle . 
The derivative terms in equation (3) have their largest value for -O' = 0 (P Q (1) =1 for all 
n when ~ 0 and cos -0 = 1). The maximum value of the terms therefore can be ob¬ 
tained from the simpler form of (3) where ^ = 0 and z =• P 1 

Hz (Z,0) - Hz ( 0 , 0 ) [l +£* (fJV £ 4 - (|j) 4 + ~] ( 5 ) 


• • 
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The most commonly used method of achieving some degree of homogeneity In coils 
is to use the spaced helrcholtz pair. Garrett has compiled a number of tables dealing 
with this method of compensation and since these tables have not been published we are 
making them part of this report. These tables are for constant current density only. 

The first table. Table 3 is used to determine the proper spacing between finite coils to 
achieve cancellation of the second derivative at the origin. The second table, Table 4 
is used to determine how far from the origin one can move before the field deviates by 
more than 0.1 percent off the field at the origin. The third table, Table 5 is particularly 
useful and can show how accurately the coils must be spaced in order to cancel the 
second derivative to the desired degree. Tho two remaining tables. Table G and Table 7 
give the resultant fourth and the sixth £ coefficients x'emaining after the proper spacing 
of the coils. The tables are used with the notation of equations (3) and (4) except that 
Garrett has normalized tbs expansion to the mean radius of tho coil rather than to the 
inner radius as In Section 13 above ( (rather than (dL ) whore ). 

To Illustrate the use of Garrett’s tables, we work out tho following examples: 


Assume we have two coils, each with an CT - 3 and a 1 and wo wish to sepa¬ 
rate them to cancel the second order derivative at the origin. 


X - 


A/i 
a*+1 


A- 


2(o-l) 
cr +1 


1.0 

- 1.0 


AX _ AX 
at+az “ atO + oT 


• • 


From Table 3, PC - - CO 








therefore 


^- 2.33040 = 2/3 4 -^ 

and the spacing should be 

'All =» . ;'3040 

a 1 

Using Table 4 we see that V- — =181 or (j+cr) ~ * 181 ^^5% T ^ s meana 

out to a distance along the axis of 72.4% of the inner radius, the field is within 0.1% of 
the -field at the origin. Using Table 5, we can examine the effect of improper spacing. 

H=H 0 (n-Q^)-Q 2 +e 4 (^) 4 4—) 

where Q is the value in table 5 and d/a Q is the displacement error in locating one coil. 
For the above example, we can find the error that will make the uncancelled 2nd deriva¬ 
tive error equal the fourth at ( — ) = 0.1 

Q -4.21 (at X “1, A =*l) 

£ 4 = -922 from table 6 
4.2t x 10" 2 = -922 x 10~ 4 

» 2.18 x 10- 3 

a 0 

D. 

When maximum homogeneity and not maximum efficiency is of the greatest impor¬ 
tance, the higher derivatives must be handled. Since the number of derivatives that can 
be cancelled depends upon the number of variables, we must now pick solutions which 
have more variables at our disposal. For Instance, if we wish to cancel both the second 
and the fourth derivative of the field we can superimpose two colls on each other each 
with two geometry variables or each with one geometiy variable and one current variable. 
To proceed, then,we would pick the length of one coil and the current of one coll and the 
variables then would be tho length of the second coil and the ratio of the two current 
densities. Garrett has calculated a set of so called sixth-order solenoids (no second or 
fourth derivatives) by maxing the following choice of variables: lie constructs an over¬ 
wound end solenoid by winding N layers of turns over a certain central region of the sol¬ 
enoid and 2N layers over a certain section of both ends (i.e. a notch on the center 
section of the o. d. of tho coll). Jlo now solves uniquely for tho length of the overall 

coll and the length of tho notch. The results of his calculation are shown In Table 8 for 
several thlcknoBsoH of coll. 


• • 



• • 
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There are other choices of two variables that ono can mako to construct sixth-order 
solenoids and Table D illustrates sevoral methods. Table 9 was constructed for the 
specific example of a ra'her long solenoid. Several general remarks can be drawn from 
the table. Probably the moat Important of which ts that since tho solcno: 


»ld la very long 







the job of compensating It for high homogeneities is much simplified, the higher derivative 
having been quite small to begin with. We notice, however, that In all but one case in 
the table, cancelling the second and fourth derivatives of the field, have increased the 
sixth derivative. The single exception 1? the current sheet of the o. d. The current sheet 
on the o.d., however, requires more power than any of the other methods of solution. 
Perhaps the most attractive method - of solution is that represented by the last column, 
that of current sheets on either end of the i.d. of the coil. This method takes little 
power and leaves the center section of the magnet open for further compensation if desired. 

Cancellation of higher derivatives than the fourth is of course possible, but is often 
of decreasing practical importance because of the unrealistic accuracy necessary In the 
cancellation of lower derivatives, if the leftover uncompensated part of say the second, la 
to be smaller than the sixth derivative or eighth derivative, etc. It probably makes sense 
to design colls to cancel the second and the fourth derivative and use a small auxiliary 
set of coils driven from the separate power source to tune out any random errors which 
creep in. Some of these random errors are of course going to be nonsymmetrical and 
will give rise to odd derivative terms. It la quite possible to construct a set of auxiliary 
compensating coils which will operate independently of each other. One coil can be 
designed to have a second derivative and no fourth and a second coll to have a fourth 
derivative and no second, etc. Each one can then be tuned independently of the other. 


We have constructed a table of <5 coefficients up to the eighth order for coils of 
varying <T and J3 with both constant current density and with current density equal 
i =* 'tf/p . The tables also include the G factor and the multiplication of the G factor 
and the & coefficients. Tbe tables are to be used with the expansion of equation (3) 
and (4) where z is referred to the inner radius Qj ; & is tabulated as D for uniform 
currents and as A for the radial current distribution. Table 10A gives the coefficients 
for the uniform current density case for several CT 3 and for /J up to 10. Table 10B 
gives the coefficients over the same range of Cf and /3 for the radial current distribu¬ 
tion case. Tables 11A, B, C, D show more detail for the constant current and the 
r adlal current cases in the small range. 


The use of the tables is probably best illustrated by the following example: Let 
Hauine wo have a coll we wish to separate for purposes of access and we want to 
' tbe resultant field in the gap and the resultant homogeneity (see Fig. 2) 











W - power 
X = space factor 
p =* resistivity 


ce ~ ^/aj 

A . 4|±4l G t (»,/*) 

^ “ 1ST Gg ( a-,y99 ) 

Figure 2 


If we now subtract from the field expansions for a coll of length 4b + JL^ . a coil of length 
we can write: 

Hgcp “ (Hot“ H Og)+(H2t“H 2g )^“)-»'(H4 t -H4 fl )^ + — 



where the primed quaiftlties refer to the error coefficients for the gap. TaVlng a spe¬ 
cific geometry, let the current be constant = D). 
let 

fit “ 4 

a* - 3 

/? g = 0.2 


What is the field In the gap and the second error coefficient, and what would be the 
second error coefficient if there were no access gap? 







The D 2 listed for CT = 3, /3 -4 In Table 10A is - .1415 x 10 - *, therefore the second 
error coefficient of the whole coll would be £z~ -1.415 x 10 ^ where 

H ~H 0 ( 1 h ^2 (57) + 

To solve for the split coO, we use (5) and the Gj. Dn columns in Table 10A 

-8-|-0 .2235 * 10" 2 4- 0.5004 x 10" 1 "*J 

H 9°P H^r(,03)[, 38 ^ 0,9 65 .O 1 ^] 

this gives a central field of 

H 9op^'°>-(^) l/Z ( 0f43 ) 

and If the second order coefficient Is written In the following form: 

H = Ho (1 + £2 (fi)** —) 

the coefficient for the above example would be 

- - 6.46 x 10~ 2 

.1 08 7 

The loss In field for this case has been 9. 5% and the second order error coefficient has 
been Increased by a factor of 4.55. If we had been trying to homogenize the field by 
separating the coils, we note that we would have separated them by too much. To ”ie 
these tables for synthesis (i. e. cancellation of derivatives) It would be necessary to use 
trial and error. In the manner of the above example. 

It is Interesting to note the oscillatory behavior of the <5 coefficients, particularly 
Lu the small yQ region. The second derivative is always negative, the fourth derivative 
goes through zero once, the sixth derivative goes through zero twice and the eighth 
derivative three times. This means that there are coils which can be built with no fourth, 
no sixth and no eighth derivatives. This may be useful In the design of Independent com¬ 
pensating coils. 








if one wishes to compensate a very long solenoid over a considerable length of the 
axis such as a multicoil structure used for plasma experiments, expansion of the deriva¬ 
tives around the origin has limited usefulness. One might wish to use the following pro¬ 
cedure instead; namely specifying that the field at any number of intervals along the axis 
be constant and then solving for ti e current distribution in the coils that would make the 
field be equal at the specified points. In the case of a multicoil structure one could ask 
that the total field in the center of each of the coils be equal to a constant and then solve 
for the current density required in each of the coils to make this true. The same proce¬ 
dure could be used with a single long solenoid where the solenoid could be broken up into 
an arbitrary number of pieces. 

The most useful method of solution of this type of problem is basically the following: 
taking any single ceil we first find the field at the center of that coil and at distances 
away from the coil represented by the center lines of each of the other coils, for a given 
current. Then by properly adding up the colls in matrix form and demanding that the 
field matrix times the current matrix be a constant we can invert the current matrix and 
solve for the individual currents required In each of the colls to achieve the desired re¬ 
sults. The method of solution and a useful notation are developed in the following example. 
Mr. R. Bradshaw of Conesco, Inc., Arlington , brought this method of solution to our 
attention. Consider the following system of coils arranged so that they have a common 
axis. They may in general be of various sizes but for illustration consider them to be 
all identical with the same separation. For ease In calculation choose an even number 
of coils so that we have a plane of symmetry passing between coils. 


Current in each coil 

Figure 3 



Our object is to find what ratio of currents is needed in the various coils so that a uni¬ 
form magnetic fie* ' Is obtained throughout the interior. If the coils are thin and close 
together, to a approximation uniformity Is obtained by requiring that the field at 

the center of e coll be the same. The field falls off rapidly on each side of the 
center of ; :. < > and we may construct a matrix the elements of which represent 

the contribution;! o. Ho colls at their respective centers. Bet Fij be the field at the 
center of coll t due to the field of coil j . For example, for 4 coils, 2 coils on each 








side of the plane of symmetry, we have 


Fn 

F 12 

IF13 

F14 

Fzt 

F22 

! F « 

*24 

f 51 

Fse 

! f 33 

F34 

F41 

F A z 

1 Ki 3 

Ft 4 


These matrix elements may be expressed in % of an Infinite solenoidal sheet of the same 
mean radius, or whatever one desires, since we need only ratios. 

F|j I j ~ 7\ where'll ls the field at the center of the coil or 


1 

F lt F i2 j Fj3 F)4 


It 


7, 

F 2 I F 22 1 F 2 3 F 2 4 

___„j_ 


h 


7 z 

F 3I F 32 j F33 F34 


h 


% 

j 

Faz j F43 F44 


I 4 


7a, 


If one now specifies «ne the Fjj matrix can be inverted and the Ij matilx 

found. 

If the current in the coils are arranged so that It ~ I 4 arid I 2 “I 3 then it is convenient 
to partition the matrix as follows: 


f.. 


fn Fie 
Fz\ Fee 


ftt- 


F13 Fj 4 
F 23 F 24 



f .. F 34 




?3 

7a 


F 11 fie ji _ 7t j and then x 2, ~ ^i| and 7z 

^21 fz2 U 7z\ where k 


O 1 
1 0 


It 

I 2 


12 


and 


fll fl2 

fei f 2z 




'Ml 


it follows that 


(f n + f T 2l<)ii - 7\ 













Since we desire the fields at the centers of each coil to be a constant 
J-\ “ 3*2 = constant = C whore C C- j | j 

i.e., + = C 

It is clear that -tfn + fntlOc and this gives us the desired ratio of currents. 


As an example consider the field on the axis due to a system of four current loops 
of radius y =» 1 Cm 



Figure 4 


Separation between centers of the coils is 1 cm. We need calculate the field of one 
loop at the center and 1, 2, and 3 cm. on the axis away from the center. For a single 
loop 

u ^ * _ y _ 

x to (x z +y2)** 


?1l- 

] 0.628 
| 0.222 

0 .222] 

0.628] 

fl2 “ 

0,056 

0.222 

0.019] 

0.056] 


] 0.056 
10.222 

0.019 I 
0.056 1 

1 ? 

O 
O - 

10.019 
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0.056] 

0 .222] 



0.628 0.222 

. 0.019 

0.056 

0.647 

0.278 

0.222 0.628 | 

0.056 

0.222 

0.270 

0.850 


We now need the inverse of ( -f*j 1 -4- 2^) 

By definition, tha inverse of a matrix A with matrix elements djj Is obtained if for 

each element (flij) -.1?__AiL ; Aji is the minor of the determinant A. 

det A 











a " a «| then A _, » 

021 a 2Zi 


For 2xZ matrices this is very simple: 


(an) '(aiz)“' 
(az.y’Cazz)-' 


det 


|_ / a 22 “ ^12^ 

tA \~a 21 an/ 


det (f t , + fizk) - (.850)(.647)-(276) 2 - .473 

ft*+W-^s 


-850 -.278 

1.797-.587 

-.278 .647 

-587 1.367 


I, 

1.797 

-587 

1 c « c I* 210 

12 

“ -.587 

1.367 

c c 

t .780 


• • 



• • 


If - c(r.2to) I 2 - c(o.7QO) 

• a 


The current In coil 1 must then be 1.55 times greater than the current in coil 2. 

Ingarden and Mlchalczyk^ have also considered the problem of homogeneities over a 

certain length of the axis rather than strictly at the origin. They have used the require- * * 

ment that the mean square deviation over a certain interval be minimized. Their notation 
is useful, but solution of the equations, which are transcendental in character, presents 
a number of problems. 

I • 


It Is often of Interest to know the field along the axis of a solenoid both Inside the 
solenoid and outside the solenoid. The following two equations can be used to find the 
field along the axis of a uniform current solenoid and a radially distributed current sole¬ 
noid (Bitter discs) respectively. 


(I) Uniform current 1c =* ^t| 

f . 7r»Xa, f (k 

5 r (k*>3M«*+(k+£) z J 


(ii) Radial current ^ 5 
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(8) 












Tabulated values of equation (7) for a wide variety of coll shapes have been published by 
D. E. Mapother and James Snyder In a University of Illinois circular #66 entitled. ’’The 
Axial Variation of the Magnetic Field In Solenoids of Finite Thickness”. There ia no 
available tabulation of equation (8j 

Equations (7) and (8) can be written in a much simpler form utilizing, the following 
angle notation: 



Figure 5 


1) uniform current density 



General expressions for the off axis fields arising from a loop of current can be 
written In terms of elliptic Integrals of the first and second kind.^ They can be written 
in the following way, where the parameters are defined in Figure 6. 


- fss [ p W sOd] 

H /° - i® [- F(k)+ (aA+ q^ z2) 3 w] 


3. 


W. R. Smythe, STATIC AMD DYNAMIC ELECTRICITY, p. 266, McGraw-Hill, New 
York, 1950 





Q - (a+/o) 2 +z 2 



F(k) 


Elliptic integral of 
the first kind 


s(k) 


EftO where E(k) = 

~ —r" 2 ~ elliptic integral 
of the second 
kind 



Using equations (11) and (12) the field at any point produced by a finite coil can be found 
by dividing up the coil Into loops and summing all the contributions. A number of labora¬ 
tories, including our own have written simple computer programs utilizing equations (11) 
and (12) to find the fields in and outside the conductor volumn. 


On-axis and off-axis fields from any solenoid with any radial current distribution 
can be calculated by means of a lengthy set of tables published by the Oak Ridge National 
Laboratory Report #ORNL 2828, entitled, "Tables for a Semi-Infinite Circular Current 
Sheet, " by N. B, Alexander and A. C. Downing. These tables are for the field around 
the end of a semi-infinite current sheet and the coil in question is built up out of these 
current sheets In the method shown in the introduction of the table. 


Attention is also called to a tabulation by E. E. Callaghan and S. H. Maslen, NASA 
D4G5 "The Magnetic Field of a Finite Solenoid". This paper plots graphs of the radial 
and axial components in and around current sheets of various lengths. While not as 
accurate as the Oak Ridge tables, it is very useful for approximate results. 


A final paper of particular relevance is that of G. R. North, "Some Parameters of 
Lumped Solenoids", (Oak Ridge ORNL-2975). He derives an expression for the field 
ripple resulting from the separation of coils in a long multicoil structure (end effects 
neglected). The equation is as follows (see Figure 7) 
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TABLE 1 

TABLE OF ERROR COEFFICIENTS 

THE £ COEFFICIENTS ARE CONSTRUCTED FROM THE 
FOLLOWING VARIABLES: 


C|“ 


c 2 = 

c 3 = 

c 4= 


1 +/J 2 

l 2 

H-/3 2 ' 

2 

cc 

a z +/3 z 

/? 2 

<r 2 +/3 2 


sin 2 £- | j 


COS 2 &- 


sin 2 ^ 0 




c 5 = im. 
C 6 = Jtn 


I <r4(/3 2 +tg 2 ) l/2 \ « / lan ('% + ° 10 /2 ) 

\ l+(/3 2 + 1) ,/2 j \ tanfe +^'/2 ) y 

/ /?+(l 2 -H)'/ 2 \ , / tan ( &w /2 ) A 

y 0+(/3 z +a z )' /z J \tan(^'i/2;y 


THE COEFFICIENTS ARE USED IN THE FOLLOWING SENSE: 

Hz = Ho (^J Pz ( u ) + <f 4^-j/°4(u) + "-j 


II “ cos# 




TABLE 1 (continued) 









TABLE I 


CASE 2 (continued) 


3/2 / 1 *3 „ ,15 2 
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TABLE 1 (continued) 


CASE 3: RADIAL CURRENT DISTRIBUTION 
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4- 4 a*ct 

, - i 

" 6 6 H b c b 
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TABLE 2 


TABLE OF LEGENDRE POLYNOMIALS 
AND THEIR DERIVATIVES (EVEN NO’S) 


p Q (u) = 1 u= COS-&- 

Pa 00= 0 

P 2 (u) = ^(3u 2 -t) 

p'z (u)=i(6u) 

/° 4 . ( u ) = -|-(55u 4 -30u 2 +3) 

/°X M = j(140u 3 ~60u) 

/°6 ( u ) = )6 ( 23 lu 6 - 315 u ^ + 105 u z - 5 ) 


/>6 (u)=i(t386u 5 -1260u s + 21C>u) 

/°8 (u) = ~g (6,435u a -12,012u 6 + 6.930u*-l,260u z +35) 


y°3 ( u ) 


128 


(5t.4S0u 7 -72.072u s +27.72Ou 3 


2,520) 






SPACEING OF HELMHOLTZ PA!R 

















AXIAL ERROR LIMITS 









TABLE 5 
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TABLE IIC 

SIXTH ERROR COEFFICIENTS 


























































